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ABSTRACT 

We have developed a model for molecular hydrogen formation under astro- 
physically relevant conditions. This model takes fully into account the presence 
of both physisorbed and chemisorbed sites on the surface, allows quantum me- 
chanical diffusion as well as thermal hopping for absorbed H-atoms, and has been 
benchmarked versus recent laboratory experiments on H2 formation on silicate 
surfaces. The results show that H2 formation on grain surface is efficient in the in- 
terstellar medium up to some 300K. At low temperatures (<100K), H2 formation 
is governed by the reaction of a physisorbed H with a chemisorbed H. At higher 
temperatures, H2 formation proceeds through reaction between two chemisorbed 
H atoms. We present simple analytical expressions for H2 formation which can 
be adopted to a wide variety of surfaces once their surfaces characteristics have 
been determined experimentally. 

Subject headings: dust, extinction - molecular hydrogen -ISM molecules - 
molecular processes 

1. Introduction 

Molecular hydrogen is the most abundant molecule in the universe and dominates 
the mass budget of gas in regions of star formation. H2 is also an important chemical 
intermediate in the formation of larger species and can be an important gas coolant in some 
conditions, particularly in the early universe. Yet, despite its importance, the H2 formation 
process is still not well understood. Obscrvationally, it has been shown that H2 can be 
efficiently formed over a wide range of temperatures ( Jura 1974; Tielens & HoUenbach 
1985a, 1985b; HoUenbach & McKee 1979). Theoretically, Gould and Salpeter (1963) 



-3- 



showed the inefficiency of H2 formation in the gas phase and concluded that recombination 
of physisorbed H on "dirty" ice surfaces is efficient between 10 and 20K. HoUenbach and 
Salpeter (1970; 1971), recognizing that this small temperature range presents a problem, 
considered reactions involving H atoms bound to dislocations and impurities with energy 
exceeding normal physisorbed energies, and obtained a recombination efficiency ~ 1 up 
to a critical temperature between 25 and 50K. Goodman (1978) calculated the quantum 
and thermal mobility of the atoms on graphite grains assuming that these atoms could 
only be physisorbed. Many studies focussed for various reasons on icy surfaces where H is 
physisorbed (Buch & Zhang 1991; Takahashi et al. 1999). However, most astrophysically 
relevant material (e.g., silicates, graphite) can bind H also in chemisorbed sites (Barlow & 
Silk, 1976, Aronowitz & Chang 1980; Leitch-Devhn & WiUiams 1984; Klose 1992; Promherz 
et al. 1993). On these surfaces, binding can occur in a physisorption layer (E~500K) at 
a distance of some Z~3A as well as in a chemisorption layer (E^IOOOOK) deeper into 
the surface Z~0.5l (Barlow & Silk, 1976; Zangwill, 1988). Recently, Katz et al. (1999) 
developed a model for H2 formation benchmarked by two sets of experiments. This model 
considers the atoms bound only in physisorbed sites and diffusing only thermally on the 
surface, colliding and recombining to form molecules. As for the ice models studied earlier, 
this allows molecule formation only below 15K for ohvine grains and 20K for carbon 
grains which contradicts the ISM observations. Perusing these models, it is clear that 
the process of the H2 formation is governed by the binding of atomic H to the surface 
and the concomitant mobility of these atoms (Leitch- Devlin & Williams 1984; Tielens & 
AUamandola 1987). Surface diffusion can occur through quantum mechanical tunneling 
(at low temperatures) and through thermal hopping (at high temperatures). For a proper 
description of molecular hydrogen formation in the ISM both types of binding diffusion 
processes have to be taken into account (Cazaux & Tielens 2002). 
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2. Model 

Recently, we have developed a model for H2 formation on grain surfaces based upon 
Langmuir kinetics where species accrete, migrate, react and the product species, possibly 
assisted by thermal energy, evaporate. This model is based on two main points: 1) The 
atomic mobilities are due to a combination of quantum tunneling and thermal diffusion 
and this is controlled by the temperature of the grain and the mass of the species. 2) 
Atoms can bind to the surface in two energetically different sites: a chcmisorption site and 
a physisorption site. These interaction energies set the migration between the different 
sites and the reactions among the species. We have compared the results of our model with 
laboratory experiments on molecular hydrogen formation on olivine surfaces (Pirronello et 
al. 1999; Katz et al. 1999) to determine the relevant parameters. Here, we extrapolate 
this model to study H2 formation under astrophysically relevant conditions (e.g., low 
accretion rates and long timescales). We consider three different populations: physisorbed 
H, chemisorbed H and physisorbed H2. The surface concentrations of these species are 
described by three rate equations. Using the surface concentration of these species, the H2 
desorption rate can be determined, which contains two contributions. First order desorption 
occurs when the grain temperature is high enough to allow evaporation of previously formed 
molecules. Second order desorption occurs when two atoms encounter and the newly 
formed molecule is directly released into the gas phase. We define Hp, Hq and as the 
physisorbed H, chemisorbed H and molecular hydrogen concentrations, respectively; their 
evaporation rate are written as IShp, Phci Ph2- The different mobilities to go from a site i 
to a site j are given by aij where i and j are either a physisorbed site (P) or a chemisorbed 
site (C), and where the barrier between two sites is assumed to be square. However, the 
shape of the barrier is found out to be unimportant (see below), n is the fraction of H2 
which stays on the surface upon formation. The rate equation can then be written as (for 
details see Cazaux & Tielens 2002), 
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Hp = F{1 -Hp- H2) - ap,Hp - 2appHp^ + a^pHcil - Hp) - (iu^Hp (1) 
He — OipcHp{l — He) — apcHpHc — o-cpHc — 2accHc^PHcHc (2) 
H2 = jJi{2appHp^ + apcHpHc + a^pHcHp + ac^Hc') - (3H2H2 (3) 

Where F is accretion rate in unit of monolayers per second (ML.s~^). The first order 
H2 desorption rate is given by: 

h = Ph,H2 (4) 

For second order desorption several processes contribute depending on whether 
physisorbed, chemisorbed atoms or both are involved: 

/C2 = (1 - iJi){apcHpHc + acpHcHp + apHp + accHc) (5) 

the desorption rate for a population i is written j3i = z^^exp [— ^] , where is the 
frequency factor of species z, perpendicular to the surface, and the desorption energy 
of species i. The total H2 formation rate, R, is the sum of the first and second order 
terms, and the H2 recombination efficiency - the fraction of the accreted hydrogen which 
desorbs as H2- is defined by e = ^ . The parameters in this model are summarized in 
Table 1. They have been discussed extensively elsewhere (see Cazaux & Tielens 2002). We 
have integrated this set of time dependent equations using a Runge-Kutta method with 
adaptative stepsize control, (Odcint, Subroutine package from Numerical Recipes 1992) 
until steady state was achieved. Under astrophysically relevant conditions, steady state 
is generally achieved. These steady state results are shown as a function of T in fig. 1 
for three different H-fluxes. All curves are characterized by a very low efficiency at low 
temperatures, a steep rise around ~5K - whose location depends on the actual H-fiux -, 
an efficiency of unity until some 25K, followed by a gradual decline to ~ 0.2 at ~ lOOK, 
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and a final drop to very low efficiency around 300K, whose location again depends on the 
actual H accretion rate. In steady state, the system of equations 1-3 can be simplified by 
setting the time derivatives equal to zero. We can now discern two different temperatures 
regimes - low temperatures (T<100K) and high temperatures (T>100K) - where different 
processes dominate H2 formation. In these limits, the set of steady state equations can be 
further simplified and this yields much insight in the numerical results. 



2.1. Physisorbed H 

At low temperatures, by far the most efficient way to form H2 is when a physisorbed 
atom reacts with a chemisorbed atom. Other routes forward H2 as well as the evaporation 
of chemisorbed atoms are negligible. Realizing that Hp is always small compare to 1 
(physisorbed H quickly migrate to a chemisorbed well), the steady state equations can be 
simplified to, 

F(l - H2) - ap,Hp - Ph^Hp = (6) 
apcHp{l - 2Hc) = (7) 
^{ap,HpHc) - Ph,H2 = (8) 

The different populations can easily be determined in this system (e.g., Hc=^, 
Hp={^^ + ^ + ^^)~^ and H2=^^Hp), the H2 desorption rate is calculated and we 
deduce the recombination efficiency: 

(s/Jijj «pc ) ^ ^ 

We can recognize three different regimes in this temperature range. At really low 
temperatures, below some lOK, the high mobility of H atoms due to tunneling permits 
the recombination of incoming H atoms. The H2 formed stays on the surface, since 
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the temperature is not high enough to allow evaporation and blocks further H-atom 

-1 

. At higher temperatures, 
(6-25K) the desorption rate depends only on the flux because all the incoming H atoms 
adsorb, recombine and desorb as H2. In this temperature regime, hydrogen accretes into 
a physisorbed site but quickly drops into a chemisorbed well. If this chemisorbed site is 
empty, the H atom will be trapped. If an H atom is already present, reaction will occur and 
the product will evaporate quickly either upon formation or thermally assisted; thus, Hc—\ 
and ejj2=l- Between ~25K up to ~100K, evaporation of physisorbed atoms competes with 
recombination, and the desorption rate decreases considerably. 



accretion. Then, the recombination efficiency is eH2= 2^ 



2.2. Chemisorbed H 

With increasing temperature, the physisorbed (at T^IOOK) and then the chemisorbed 
atoms(at T^^SOOK) will start to evaporate. The most efficient reaction to form H2 is then 
the "coUision" of two chemisorbed atoms. The system of equations reduces to, 

F - apcHp - pHpHp = (10) 

apcHp - 2a,,Hl - (3hcHc = (11) 
pL{a,,Hc^) - Ph,H2 = (12) 

and the H2 recombination efficiency can be written as. 



1 + ^) \ (13) 



with ^ the correction factor at high temperatures 



'2F OtpcOicc 



(14) 



which reflects the evaporation of chemisorbed H. This factor is flux dependent. 
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2.3. General expression for the recombination efficiency 

The above expression for the H2 recombination efficiency can be combined into one 
general expression valid at any temperatures: 



The mobilities, apc and acc-i ^-re due to a combination of quantum mechanical tunneling 
and thermal hopping but considering the temperatures at which these parameters play a 
role, these mobilities are dominated by thermal hopping. Therefore, the H2 recombination 
efficiency is independent of the width of the considered barriers, and we can approximate ^ 
and by the following expressions: 



1 + 



phc exp 



l.bEHr 



kT 



2F 



-1 



(16) 



a 



pc 



Ehc - Es 
Erp - Es 



exp 



Es 
kT 



(17) 



The parameters in these expressions have been determined from experimental data 
(Table 1; Cazaux & Tielens 2002). The H2 recombination efficiency for three different flux 
is reported figure 1. In astrophysical environments, the recombination rate is written: 



= -jnHVHnd<JdeH2SH{T) 



(18) 



where ur and vh are the number density and the thermal velocity of H atoms in the gas 
phase, ridOd is the total cross section of the interstellar grains and Sh{T) is the sticking 
coefficient of the H atoms which can depend on temperature. 
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3. Discussion 

Our study reveals the presence of two distinct regimes of H2 formation, which reflect 
directly the presence of two types of atomic H binding sites. At low temperature (T<100K), 
H2 formation involves migration of physisorbed H atoms. At higher temperatures 
(T>100K), H2 formation results from chemisorbed H recombination. The presence of 
these two types of binding sites allow H2 formation to proceed relatively efficiently even 
at elevated temperatures. The study of Hollenbach and Salpeter (1971) focused on icy 
surfaces on which H can only physisorb. As a result, H2 formation ceased at temperatures 
in excess of ~20K. Recognizing this problem, Hollenbach and Salpeter involved the presence 
of enhanced binding sites on the ice with ill-determined parameters. These sites allowed H2 
formation to proceed up to some 75K. Since their study, it has become abundantly clear 
that interstellar grains are not covered by ice in the diffuse interstellar medium (Whittet et 
al. 1983, 1988). Sihcate and graphitic surfaces are now widely accepted as astrophysically 
relevant grain surfaces (Mathis 1991) and those surfaces intrinsically possess enhanced 
binding sites; e.g., chemisorbed sites. The parameters of these chemisorbed sites have 
not yet been well determined because experiments have focused on low temperature H2 
formation (Pirronello et al. 1997; Katz et al. 1999). The values adopted in this study 
are however quite reasonable and illustrate the efficiency of H2 formation at elevated 
temperatures well. When future experiments determine the values of the parameters 
involved {Ehc, Ehp, jj), the results can be directly adjusted. Inside dense clouds, interstellar 
grains are covered by ice. Of course, in such environments, almost all hydrogen is already 
in molecular form and H2 formation is perhaps only of academic interest. Nevertheless, 
we note that in such environments, molecular hydrogen formation may proceed mainly 
through H-abstraction from molecules such as H2S and N2H2 (Tielens and Hagen 1982). In 
a sense, these species act as "chemisorption" sites for hydrogen. Migrating atomic H may 
tunnel through the reaction barriers involved and form H2. Eventually, these ice covered 
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grains are transported back into the diffuse ISM, wfien tfie molecular cloud is disrupted. 
Photodesorption and sputtering in strong shocks quickly remove their ice on a timescale of 
some 10^ yr (Tielens and Hagen 1982, Jones et al. 1994, Draine and Salpeter 1979, Barlow 
1978). At that point, molecular hydrogen formation is again governed by the properties of 
bare grain surfaces. Similarly, any (thin) layer accreted in the diffuse interstellar medium 
will be quickly sputtered in even a modest velocity (~30km s~^) shock (Jones et al. 1994). 
Finally, the formation efficiency of molecular hydrogen will also depend on the sticking 
coefficient of H atoms colliding with the grain. In our model and the formulae derived, the 
sticking coefficient is subsumed in the incident flux, F. Astrophysical studies of sticking 
of H on grain surfaces have concentrated on physisorbed interactions and the sticking 
coefficient is ~1 at low temperatures and decreases with increasing temperature to about 
0.3 at T=300K ( HoUenbach & Salpeter 1970; HoUenbach & McKee 1979; Burke & 
HoUenbach 1983; Leitch-Devlin & Williams 1985). However, if the interaction occurs 
through much stronger chemisorption, then the sticking coefficient might be large even at 
high temperatures (Tielens & AUamandola 1987; Duley & WiUiams 1984). 



3.1. Summciry and conclusions 

Recently, we have modelled molecular hydrogen formation on grain surfaces. This 
model consider hydrogen atoms bound to the surface at two energy levels (i.e. chemisorption 
and physisorption) . The H mobility from one site to another is a combination of tunelling 
effect and thermal diffusion. This model has been experimentally benchmarked (Pirronello 
et al. 1997a, 1997b, 1999) and the relevant surface characteristics have been determined. 
These characteristics allow us to extend our model for H2 formation under astrophysically 
relevant conditions. The results show efficient H2 formation from ~6K to ~300K. The 
different processes involved in H2 formation at different temperatures has been discussed. 
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Until about lOOK, H2 forms by recombination of a physisorbed H with a chemisorbed H 
and is highly efficient. At higher temperatures, when physisorbed atoms evaporate quickly, 
the recombination of two chemisorbed atoms is required to form H2. H2 formation is then 
less efficient, ei72~0.2. The parameters involved in H chemisorption and H2 formation at 
high temperatures are presently not well known. The adopted values are very reasonable 
and the gross characteristics - H2 formation at high temperatures - are undoubtedly 
correct. Nevertheless, future experiments are very important to determine the maximum 
temperature to which H2 formation in the ISM can occur. 
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Table 1: Model parameters for silicate surface. 
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320 


0.005 


200 


600 


10000 


3 10^2 


1.3 10^3 



for more details about the determination and calculation of these parameters, see Cazaux & Tielens 2002 

Note. — Eh2, Ehp and Erc are the desorption energies of H2, physisorbed H (Hp) and chemisorbed H 
(He), and Es is the energy of the saddle point between two physisorbed sites, fi is the fraction of the newly 
formed H2 which stays on the surface and ur^ and vrc ai'e the vibrational frequencies of H2 and H in their 
surface sites. 
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Fig. 1. — H2 recombination efficiency for tfiree different fiux: 10 monolayers per second 
(solid), 10~^° monolayers per second(dasli) and 10~^ monolayers per second (dot-dash). 



